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Abstract Sympatric reptiles are the ideal system for investigating temperature-driven coexistence. Understanding 
thermally physiological responses of sympatric lizards is necessary to reveal the physiological mechanisms that underpin 
the sympatric occurrence of reptiles. In this study, we used three lizard species, Eremias argus, E. multiocellata, and 
Phrynocephalus przewalskii, which are sympatric in the Inner Mongolia desert steppe, as a study system. By comparing 
their resting metabolic rates (RMR) and locomotion at different body temperatures, we aimed to better understand 
their physiological responses to thermal environments, which may explain the sympatric occurrence of these lizards. 
Our results showed that Е. argus had significantly higher RMR and sprint speed than Е. multiocellata, and higher 
RMR than P. przewalskii. In addition, the optimal temperature that maximized metabolic rates and locomotion for E. 
argus and Е. multiocellata was 36°C, whereas for Р przewalskii it was 39°C. Our study revealed the physiological 
responses to temperatures that justify the sympatric occurrence of these lizards with different thermal and microhabitat 
preferences and active body temperatures. Eremias argus and Е. multiocellata, which have lower body temperatures 
than P. przewalskii, depend on higher RMR and locomotion to compensate for their lower body temperatures in 
field conditions. Our study also highlights the importance of using an integrative approach, combining behavior and 
physiology, to explore the basis of sympatric occurrence in ectothermic species. 
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1. Introduction food preferences (Hardin, 1960; Shurin et al., 2004). 


Pacala and Roughgarden (1985), for example, reported 


Understanding the mechanisms that allow species 
coexistence is one of the core issues in community 
ecology. Niche differentiation is considered to be the 
basis of coexistence among sympatric species. Niche 
differentiation is a process by which competing species 
utilize environmental resources differently, and it can 
include aspects such as activity period, use of space, and 
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the existence of several anole lizards that shared food 
resources (i.e., insects) on the Caribbean islands but 
occupied different microhabitats such as leaf litter floor 
and branches to avoid any possible competition for 
microhabitats or food resources. Lophophorus sclateri and 
Ithaginis cruentus, two pheasant species, were observed 
to share high-altitude habitats in Gaoligong Mountain but 
to feed on different plants or different parts of the same 
plants (Luo et al., 2016). 

In ectotherms, performances of specific activities such 
as locomotion, immunity, growth, and even reproduction 
are greatly affected by body temperatures (Angilletta, 
2009; Angilletta et al., 2002). Within thermal tolerance 
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range, the physiological performance of ectotherms is 
known to enhance as body temperature increases until the 
optimal temperature is reached; after that, performance 
rapidly decreases if body temperature continues to 
increase (Huey and Kingsolver, 1989). Because of its 
effect on body temperature and thus on function and 
performance, environmental temperature is one of the 
most important ecological factors for ectotherm animals 
(Angilletta, 2009) and, therefore, is also considered 
an important resource (e.g., Li et al., 2017). Among 
sympatric ectotherms, different temperature preferences 
may result in different body temperatures, which in turn 
may result in segregation of microhabitat utilization and 
allow sympatric occurrence (Adolph, 1990; Hertz, 1992; 
Martinvallejo et al., 1995; Wilkinson and Grover, 1996). 
Understanding how variation in thermal environment 
affects sympatric species has become a hot topic in 
animal ecology (e.g., Li et al., 2017; Osojnik et al., 2013; 
Rusch et al., 2018; Zagar et al., 2015). 

As ectotherms, sympatric lizards constitute ideal 
systems for investigating temperature-driven niche 
differentiation and coexistence (Pianka, 1986). Lizards 
can regulate their body temperatures in a narrow range 
mainly through behaviors to facilitate physiological 
functions (Adolph, 1990; Angilletta, 2009). On the 
Mongolian plateau, three ground-dwelling sympatric 
lizard species coexist in arid, semi-arid, or grass 
lands area: Eremias argus, E. multiocellata, and 
Phrynocephalus przewalskii (Zhao et al., 1999). These 
sympatric species have been demonstrated to occupy 
different microhabitats (e.g., vegetation coverage) and to 
have significantly different thermal preferences and active 
body temperatures. Phrynocephalus przewalskii mainly 
selects open and warm microhabitats and has a higher 
preferred body temperature range (T.a: 33.9-39.2°C ) than 
Е. argus (Та: 32.8-37.5°C) and Е. multiocellata (Та: 
33.4—36.8°С), which occupy shaded and cool habitats 
(Table S1) (Li et al., 2017). Accordingly, the active body 
temperatures are significantly higher in P. przewalskii 
than іп Е. argus and Е. multiocellata (Table S1) (Li et al., 
2017). These species’ demands for thermal environmental 
resources were speculated to be met because of 
microhabitat differentiation, which may be an important 
basis for their sympatric occurrence (Li еѓ al., 2017). 
Nevertheless, as ectotherms, their body temperatures 
should be effective on fitness-related physiological 
functions before they affect their sympatric occurrences 
(Angilletta, 2009; Angilletta et al., 2002; Hochachka and 
Somero, 2002). Our current knowledge on these species’ 
differences regarding thermal preference or active body 
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temperatures is insufficient to explain the effect of their 
thermal niche-partitioning on promoting their sympatric 
occurrence. Therefore, in order to reveal the physiological 
basis that underpins their temperature-driven sympatric 
occurrence, it is critical to investigate the effects of body 
temperature on important physiological traits of Е. argus, 
Е. multiocellata, and P. przewalskii. 

Metabolism is one of the most important physiological 
processes because it determines an organism’s demands 
from the environment and energy allocation among 
functions (Brown et al., 2004; McNab, 2002). Lizards’ 
metabolic rates have been demonstrated to be significantly 
sensitive to body temperatures; these, in turn, require 
the animals’ fast response to any environment thermal 
variation (e.g., Ma et al., 2018a; Ma et al., 20186; Sun 
et al., 2018). For example, the metabolic rates were 
enhanced as body temperature increasing, and could 
respond to thermal variation by acclimation in lizard 
(e.g., Sun et al., 2018). Furthermore, the segregation of 
metabolic rates at the same temperatures may help niche 
differentiation (Žagar et al., 2018; Žagar et al., 2015). In 
addition, locomotion is thermal dependent and critical 
to escape danger, forage, and choose mating partners, 
and thus affects survival and reproduction in ectotherms 
(Robson and Miles, 2010; Shine, 2003; Shu et al., 2010; 
Wilson, 2001). 

Here, with sympatric Е. argus, Е. multiocellata and 
P. przewalskii from Shierliancheng of Inner Mongolia 
area as a study system, we measured the metabolic 
rates and locomotion of each species at a range of body 
temperatures from 15—42°C in order to determine inter- 
specific differences in physiology and fitness-related 
performances. We predicted that: 1) Е. argus and Е. 
multiocellata, which occupied shaded microhabitats and 
have lower active body temperatures, would have higher 
metabolic rates and sprint speeds but lower optimal 
temperatures for sprint speed than P. przewalskii at a same 
temperature within a moderate range before the optimal 
temperature is reached; 2) combined with temperature- 
related curves of metabolic rates and sprint speed, the 
three sympatric species would perform better when at 
their own active body temperatures under preferred 
microhabitat conditions. 


2. Materials and Methods 


2.1. Lizard collection and husbandry The ground- 
dwelling lizards E. argus, E. multiocellata, and P. 
przewalskii were collected from late June to early July 
at the Shierliancheng field station, Institute of Grassland 
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Research of the Chinese Academy of Agricultural 
Sciences (111°09' Е, 40°21’ М, elevation 1010-1021 
т), where the ground-dwelling lizards (1.е., Е. argus, Е. 
multiocellata, and Р. przewalskii) sympatry (11 её al., 
2017). During collection, gravid female lizards were 
determined by palpate and released. Male and nongravid 
female adults of each species were transferred to our 
laboratory in Beijing, where they were weighed (to the 
nearest 0.001 g) and measured (to the nearest 1 mm). 
Every five or six lizards of each species were kept in 
a terrarium (550 х 400 х 350 mm’, length х width х 
height). The terraria were placed in a temperature- 
controlled room at 20°C with photoperiod cycles of 10 
h dark: 14 h light (6:00-20:00). A heating bulb (50 w) 
was placed above one end of each terrarium to provide a 
temperature gradient inside the terrarium from 20—50°C 
within each photoperiod. Food (crickets dusted with 
vitamins) and water were provided ad libitum. 


2.2. Resting metabolic rate After seven days of 
husbandry, the resting metabolic rates (RMR) of 14 Е. 
argus (7 3 and 7 ©), 13 Е. multiocellata (5 3 and 8 $), 
and 14 Р przewalskii (7 3 and 7 $) (total n = 41) were 
measured at ten test temperatures (15, 18, 21, 24, 27, 
30, 33, 36, 39, and 42°C), with one temperature every 
other day in a random sequence. Lizards were fasted 
for at least 12 hours before each measurement. Initially, 
each lizard was placed in an incubator (Sanyo, Japan) 
at the selected test temperature for 90 min; then, the 
lizard was placed in a respirometry chamber within the 
incubator and the measurement was performed. RMR was 
measured using a closed-flow respirometry system with 
a volume of 281.4 ml (Stable Systems International Inc. 
Las Vegas, NV, USA) and estimated via CO, production 
rate following a previously established method (Sun et 
al., 2018). In brief, the system contained the following 
modules: Universal interface II (UI2), Subsampler TR-S 
(SS3), Mass flow control unit (MFC-2), Oxygen analyzer 
(FC-10a), and Carbon dioxide analyzer (Ca-10a). At the 
beginning of each measurement, the system was opened 
for two to three minutes so that air came through an 
entrance tube with a flow rate of 300 ml/min to make 
the baseline stable. Afterward, the measurement system 
was transferred to a closed-circuit respirometry, and 
then carbon dioxide production rates (Исо) in closed- 
circuit were continuously recorded for 10 min. The 
entire environment was dark (no light exposure) during 
measurements, and all measurements were conducted 
from 10:00-18:00 to minimize the effects of circadian 
rhythms. The metabolic rates were calculated as the CO, 
production per gram”” of body mass рег hour (ml/g’”*/ 


hr) following the ‘metabolic theory of ecology’ (Brown et 
al., 2004), with the equation of metabolic rates = Исо X 
volume/body mass, where Vo, is the CO, production rate 
in percentage (%/hr) in the closed circuit with the volume 
of 281.4 ml (Sun et al., 2018). After each measurement, 
lizards were transferred back to the terraria until the next 
measurement initiated. After around 20 days of RMR 
tests, the locomotion tests were performed. 


2.3. Locomotion The locomotion of 37 lizards was 
determined using a 1 000 mm custom-made race track 
at ten test temperatures randomly (15, 18, 21, 24, 27, 
30, 33, 36, 39 and 42°C) after the RMR measurements. 
Locomotion measurements were not performed in four 
E. argus as they were out of conditions (i.e., inactive 
or unwilling to move under stimulation) one day 
after the RMR experiments. Before the locomotion 
measurement, lizards were placed in an incubator at each 
test temperature for 90 min. The sprint speed was then 
determined in a 1 000 mm race track with photoelectric 
timers every 200 mm. For measurements, the lizard was 
placed in one end of the track and then stimulated on the 
tail with a paintbrush to run along the race track. The time 
spent by the lizard to run over every 200 mm interval 
was recorded using the photoelectric timer. Tests were 
conducted from 10:00-14:00. Each individual was tested 
twice with an interval of one hour between tests; the 
fastest speed of the two tests (i.e., 10 speed records) was 
recorded as the sprint speed for each lizard. Tests were 
conducted every other day. 


2.4. Statistical analysis Before statistical analysis, 
Kolmogorov-Smirnov test and Levene’s test were 
conducted to detect data normality and variance 
homogeneity. Repeated-measures ANOVAs were 
conducted to determine species differences in RMR 
and locomotion, with species as a main factor and test 
temperatures as a repeated factor. When interaction 
between species and test temperatures was detected, 
a further comparison was conducted to analyze the 
differences among species at relevant test temperature 
ranges. 


3. Results 


3.1 Resting metabolic rates The RMRs were 
significantly different among the three species. Eremias 
argus had a higher RMR than Е. multiocellata and P. 
przewalskii (Е. argus* > Е. multiocellata®> Р. przewalskii’; 
Е, з= 3.65, Р = 0.036) (Figure 1). RMR increased 
with test temperature until the optimal temperature was 
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reached; then, RMR decreased as temperature continued 
to increase (repeated factor ‘test temperature’: Fy 34) = 
97.956, P < 0.0001). The optimal temperature of RMR 
for E. argus and E. multiocellata was 36°C, whereas 
for P. przewalskii it was 39°C. The effect of the test 
temperatures on RMR was species-dependent. At low 
temperatures (from 15—36°C), Е. argus had a significantly 
higher RMR than Р przewalskii, with Е. multiocellata in 
between (Е. argus* > Е. multiocellata” > Р. przewalskii”; 
Е, зз = 8.251, Р = 0.001). At high temperatures (39-42°C), 
the RMR of P. przewalskii was significantly higher than 
those of E. argus and Е. multiocellata (Р. przewalskii 
>E. argus” > Е. multiocellata’; Рз = 18.023, Р < 0.0001) 
(Figure 1). 


3.2. Locomotion Sprint speed of E. argus was 
significantly higher than that of Е. multiocellata, whereas 
sprint speed of P. przewalskii was not significantly 
different from those of either species (Е. argus* > P. 
przewalskii® > Е. multiocellata’; Е, 4 = 5.376, Р = 0.009) 
(Figure 2). Test temperature had a significant effect on 
sprint speed (Fo 39, = 61.687, Р < 0.0001). Sprint speed 
increased with temperature until the optimal temperature, 
after which it decreased as temperature continued to 
increase. The optimal temperature of sprint speed for 
E. argus and E. multiocellata was 36°C, whereas for P. 
przewalskii it was 39°С. At low temperatures (15—36°С), 
the sprint speed of Е. argus was significantly higher than 
that of Е. multiocellata; the sprint speed of P. przewalskii 
was not significantly different from those of either species 
(Е. argus* > Р. przewalskii” > Е. multiocellata’; Е, за 
= 4.439, Р = 0.019). At high temperatures (39-42°C), 
the sprint speeds of Р. przewalskii and Е. argus were 
similar and both higher than that of Е. multiocellata (Р. 
przewalskiť > Е. argus* > Е. multiocellata’; Е, «= 6.792, 
Р = 0.003) (Figure 2). 


4. Discussion 


On the basis of previously known thermal biology traits 
such as active body temperatures, and thermal preference 
of the sympatric lizards Е. argus, Е. multiocellata, and 
P. przewalskii (Tables S1) (Li et al., 2017), in the present 
study we determined the interspecific differences in RMR 
and locomotion at different body temperatures from 15 
to 42°C. We found that Е. argus have significantly higher 
RMR and locomotor performance, and Е. multiocellata 
has significant higher RMR, when compared to P. 
przewalskii, especially at low body temperatures 
from 15 to 36°C, but lower optimal temperatures for 
RMR and locomotion, which indicates the existence 
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Figure 1 Resting metabolic rate (RMR) of Eremias argus, E. 
multiocellata, and Phrynocephalus przewalskii at 15, 18, 21, 24, 27, 
30, 33, 36, 39, and 42°C. RMR was expressed as CO, production 
рег g°” body mass per hour (ml/g”’*/hr). Red circles, blue triangles, 
and green rectangles indicate the RMR of E. argus, E. multiocellata, 
and P. przewalskii, respectively. The optimal temperatures for E. 
argus, Е. multiocellata, and P. przewalskii were 36°C, 36°C, and 
39°C, respectively. Data are shown as mean + SE. The RMRs 
(expressed as ml/g/hr) are also provided in Table $2, for the 
convenience of interspecific comparison in ‘Meta-Analysis’. 
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Figure 2 Locomotion of Eremias argus, Е. multiocellata, and 
Phrynocephalus przewalskii at 15, 18, 21, 24, 27, 30, 33, 36, 39, and 
42°C . Locomotion was expressed as sprint speed (m/s). Red circles, 
blue triangles, and green rectangles indicate the sprint speed of Е. 
argus, Е. multiocellata, and P. przewalskii, respectively. The optimal 
temperatures for E. argus, Е. multiocellata, and P. przewalskii were 
36°C , 36°C , and 39°C , respectively. Data are shown as mean + SE. 


of a physiological adaptation to the decrease in body 
temperature that also influences their performances. 
Therefore, we found that physiology and performance are 
both fine-tuned to the thermal preferences of the studied 
lizard species. 

The interspecific differences in metabolic rates 
among Е. argus, Е. multiocellata, and P. przewalskii are 
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consistent with the differences that were observed in 
hatchlings incubated under fluctuant temperatures (Ma 
et al., 2018a, 2018b). The metabolic rates of hatchlings 
of Е. argus and Е. multiocellata were higher than in those 
of P. przewalskii and enhanced with the increase in the 
test temperatures; this may indicate that the interspecific 
differences are fixed at the different life cycle stages of 
these species. As the fundamental physiological process, 
metabolic rates may reflect the organism’s physiological 
response to the environment (Brown et al., 2004; 
McNab, 2002). However, unlike endotherms—which 
have thermal neutral zones within which organismal 
metabolic rates are thermally insensitive—, ectothermic 
metabolic rates are significantly thermally dependent, 
with a slow increasing rate at low temperatures and a 
rapid increase at high temperatures, followed by a steep 
drop after optimal temperature is reached (Gillooly et 
al., 2001; White, 2011). As the fundamental ‘pacemaker’ 
of biological rate, metabolic rate is related to lizard 
species’ body temperature, which is in turn affected 
by thermal environments (Brown et al., 2004; Glazier, 
2015). Eremias argus and E. multiocellata select shaded 
microhabitats and thus have lower body temperatures 
than Р. przewalskii (Table S1) (Li et al., 2017). Having 
high metabolic rates at low temperatures (15—-36°C, 
Figure 1) may create advantages for Е. argus апа Е. 
multiocellata at low temperatures by allowing the 
allocation of more metabolic energy for the performance 
of activities such as escaping or foraging (e.g., Sun et al., 
2018; White et al., 2012; White and Kearney, 2013). In 
contrast, as a higher RMR may enable a higher metabolic 
energetic production, a higher RMR at high temperatures 
in P. przewalskii may be responsible for improving this 
species’ performance, accompanied by higher thermal 
preferences and active body temperatures. Alternatively, 
a higher RMR at high temperatures may induce more 
energetic allocation for maintenance, which might be a 
cost (Sokolova et al., 2012). 

The existence of interspecific differences in locomotion 
across body temperatures has been demonstrated in 
numerous species of lizards (Chen et al., 2003; Du et 
al., 2000; Ji, 1995; Ji et al., 1996; Sun et al., 2014; Xu, 
2001; Zhang and Ji, 2004). Locomotion could reflect the 
species’ ability to escape, forage, and even reproduce 
(Husak and Fox, 2006; Shu et al., 2010). The higher 
sprint speed in Е. argus than in Р. przewalskii and Е. 
multiocellata may result in selective advantages in 
escaping, foraging, and even reproduction, especially 
at low temperatures (Bergmann and Irschick, 2010; Xu, 
2001; Zamora-Camacho et al., 2014). 


As they are desert species, the optimal temperatures 
for sprint speed in Е. argus (36°C), Е. multiocellata 
(36°C), and P. przewalskii (39°C) are higher than those 
reported for skinks (32—34°C ) (Du et al., 2000; Ji, 1995; 
Xu, 2001) and grass lizards (28—-34°C ) (Chen et al., 
2003; Ji et al., 1996; Zhang and Ji, 2004). In addition, the 
high optimal temperature for P. przewalskii found in this 
study is consistent with this species’ higher active body 
temperatures, thermal preference, and thermal tolerance 
(Li et al., 2017) if compared to the sympatric Е. argus and 
E. multiocellata. The differences in optimal temperatures 
may be an evolutionary consequence of adaptation to 
different climatic environments (Zhang and Ji, 2004). 
Lizards in open and warm environments (e.g., deserts) 
tend to have high body temperatures, and their optimal 
temperatures for functions are normally positively 
correlated with body temperatures (e.g., Ji, 1995; Zhang 
and Ji, 2004). Potentially, active body temperatures and 
optimal temperatures for function may be both affected 
by habitat thermal environments through natural selection 
or acclimation (Gilbert and Miles, 2017; Sinclair её al., 
2016). Within a species’ distribution area, other climatic 
factors (1.е., humidity) may also directly or indirectly 
induce variation in thermal traits in lizards, including 
optimal temperatures of thermal performance curve or 
critical temperatures (Sinclair et al., 2016; Sunday et 
al., 2012). Nonetheless, in the present study, these three 
sympatric species lived in adjacent microhabitats with 
very similar precipitation and evaporation indices (Wang 
et al., 2016); thus, we propose that the divergences in 
optimal temperatures of metabolic rates and locomotion 
among species were driven by thermal environments. 
However, future studies on the effect of multiple factors 
on thermal traits in sympatric ectotherms systems would 
be very important and necessary. 

In general, the sympatric species Е. argus, Е. 
multiocellata, and P. przewalskii have different preferred 
microhabitats, resulting in interspecies variation of 
thermal preferences and active body temperatures (Table 
51) (Li et al., 2017). Eremias argus and Е. multiocellata 
resort to high physiological processes within a low 
temperatures range (1.е., 15-36°С) to compensate 
for their lower body temperatures in field conditions. 
Alternatively, Р. przewalskii occupies open and warm 
microhabitats, and it has higher body temperatures. Lower 
metabolic rates and sprint speed at low body temperatures 
may be an adaptive strategy for Р przewalskii; its low 
energetic production at low body temperatures results in 
a lower energetic cost—indicated by low sprint speed—, 
which іп turn reduces energy expenditure (Brown et al., 
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2004; Young et al., 2011). Similarly, within a certain 
geographical scale, populations or species that occupy 
warm habitats (e.g., tropical region) also have higher 
body temperatures, and therefore tend to have lower 
performance (e.g., swimming speed) and metabolism 
rates when compared to species from cold environments 
at even temperatures, as ‘Metabolic Cold Adaptation’ 
predicts (e.g., White et al., 2012). Given that, the 
combination of different active body temperatures and 
physiological trait responses is an effective solution 
that allows the sympatric occurrence of the three desert 
lizards, by enabling their functions within an active body 
temperature range. Our study highlights the importance 
of integrative investigations on temperature-driven 
sympatric occurrence at a physiological level, based on 
thermal biological traits and active body temperatures. 
Future studies should focus on the dynamic of the effects 
of body temperatures on physiological responses in field 
to reveal the modification of physiological and body 
temperatures to thermal variation. Especially in the 
context of climate warming, studies integrating body 
temperatures alterations and thermal performance curves 
may provide insight into the responses and the evaluation 
of the vulnerabilities of sympatric species (Sinclair et al., 
2016). 
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Appendix 


Table S1 Average active body temperatures (7,), active body temperatures at open microhabitat (T, open), active body temperatures at 
shaded microhabitat (T, shade), and thermal preference (T) of Phrynocephalus przewalskii, Eremias argus, and Е. multiocellata. Lizards 
were collected at the Shierliancheng field station, Institute of Grassland Research of the Chinese Academy of Agricultural Sciences (111°09' 
+ SE and follow Li et al. (2017). 


Е, 40°21'N, elevation 1 010—1 021 т). Data are expressed as mean 4 


N P. przewalskii N E. argus N Е. multiocellata 

T, (CC) 129 38.1 = 0.2 84 37.5 = 0.2 116 36.6 = 0.2 
Т, ореп (°С) 43 40.19 + 0.88 29 39.07 = 0.67 39 38.33 = 0.14 
T, shade (°C) 40 35.83 + 1.78 27 35.70 = 1.69 32 34.81 +0.19 

36.6 = 0.2 35.3 = 0.2 35.2 + 0.2 
Та °С) 17 21 21 

33.9—39.2 32.8—37.5 33.4-36.8 
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